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Abstract. The aim of this paper is to estabhsh time decay properties and disper- 
^ . sive estimates for strictly hyperbohc equations with homogeneous symbols and with 

' time-dependent coefScients whose derivatives belong to L^(M). For this purpose, 

the method of asymptotic integration is developed for such equations and repre- 
sentation formulae for solutions are obtained. These formulae are analysed further 
to obtain time decay of L^-L'^ norms of propagators for the corresponding Cauchy 
problems. It turns out that the decay rates can be expressed in terms of certain 
geometric indices of the limiting equation and we carry out the thorough analysis 
of this relation. This provides a comprehensive view on asymptotic properties of 
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(~| , solutions to time-perturbations of hyperbolic equations with constant coefficients. 

Moreover, we also obtain the time decay rate of the L^-L'^ estimates for equations 
of these kinds, so the time well-posedness of the corresponding nonlinear equations 
with additional semilinearity can be treated by standard Strichartz estimates. 



1. Introduction 



(N 
> 
00 

I This paper is devoted to several aspects of strictly hyperbolic equations of higher 

^ ■ orders or of strictly hyperbolic systems with time-dependent coefficients. In particu- 

lar, we will investigate the following topics: 
^ ■ • representation of solutions of equations of higher order; 

! • dispersive estimates for solutions. 

Equations of orders larger than two appear often in the analysis of large first order 
systems and in the analysis of coupled equations of higher orders. In the present paper 
^ ' we will restrict our attention to the investigation of equations with homogeneous 

■ symbols (at the same time making a suitable preparation for the further development 

of this topic for equations with low order terms). In fact, we will concentrate on 
scalar equations of some order m G N keeping in mind that in the case of a system 
its dispersion relation (the determinant) will be of such form, so the information on 
solution to the Cauchy problem for the dispersion relation will imply the information 
on solutions to the Cauchy problem of the original system. On one hand, we will 
introduce several techniques allowing to deal with equations of higher orders. On the 
other hand, already for the second order equations the new method that we propose 
in this paper will yield certain improvements and extensions of known results. In 
particular, we will improve the result on the decay rates in the dispersive estimates 
already for the time-dependent wave equation, as well as the time decay rate for 
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the standard Kirchhoff equation, thus also improving the corresponding Strichartz 
estimates. It will also allow the inclusion of mixed terms in second order equations (a 
question which is known to be very delicate if we want to treat problems outside the 
perturbation framework). We will allow time-dependent coefficients and will assume 
that their derivatives are in L^(R). It is known that this property is satisfied in 
many situations, for example in applications to Kirchhoff equations and systems, 
etc. This will also allow us to obtain a comprehensive view on time-perturbations of 
equations with constant coefficients, in which case the assumption of the integrability 
of derivatives of coefficients is quite natural. For the purposes of this paper, we will 
develop the asymptotic integration method for hyperbolic partial differential equations 
with time-dependent coefficients. While this method is relatively well-known in the 
theory of ordinary differential equations (see e.g. Hartman [7]), its use in partial 
differential equations appears to be new. 

We note that equations with constant coefficients have been thoroughly studied by 
Sugimoto in a series of papers [2ll [25l |26] who described several interesting geometric 
quantities responsible for the rate of the time decay of U'-L'^ norms of their propa- 
gators. In particular, one has to look at the level sets of the characteristic roots of 
the symbol and at the orders with which tangent lines touch these sets. These orders 
become responsible for the time decay rate in the corresponding dispersion estimates 
and for indices of the subsequent Strichartz estimates. In fact, the appearing indices 
are related to the oscillation indices of integral kernels of the propagators, viewed as 
oscillatory integrals, and their classification is well studied in the singularity theory 
(e.g. [U). 

The case of strictly hyperbolic equations with constant coefficients with lower order 
terms has been thoroughly investigated in [22]. In particular, properties of charac- 
teristic roots are crucial in determining exact decay rates and the complete analysis 
is quite lengthy and involved. For example, in the case of equations of dissipative 
types analysed in [21] the decay rate is determined by properties of characteristics for 
small frequencies. A general analysis of this type is necessary for application to large 
systems, such as Grad systems in gas dynamics, or to Fokker-Planck equations, in 
which case the Galerkin approximation produces a sequence of scalar equations with 
orders going to infinity, see e.g. [19j. Applications to such problems give a strong 
additional motivation to the investigation of equations of higher orders of the type 
of those treated in this paper. 

To become more precise, we consider the Cauchy problem for an irfi^ order strictly 
hyperbolic equation with time-dependent coefficients, for function u = u{t,x): 

(1.1) L{t,Dt,D,)u = D^u+ ci,j{t)D^^Diu = 0, t^O, 

\i/\-\-j='m 

J < m — 1 

with the initial condition 

(1.2) D^u{0,x) = fk{x) e C^iR''), A; = 0,1,... ,m-l, x G M", 

where A = -idt and = {-id^^Y^ ■ ■ ■ {-id^J"", i = for z/ = (z/i, . . . , z/„). 

Denoting by Lipioj.(M) the space of all functions which are locally Lipschitz on R, we 



DISPERSIVE ESTIMATES FOR HYPERBOLIC EQUATIONS 



3 



assume that each a^jit) belongs to LipiQ^(M) and satisfies 
(1.3) G L\W) with \v\+j = m. 

Moreover, following the standard definition of equations of the regularly hyperbolic 
type (e.g. Mizohata [U]), we will assume that the symbol of the differential operator 
L(t, Dt, D^) has real and distinct roots ipiit^S,), ■ ■ ■ , ^Pm(t] ^ 7^ O5 that 

(1-4) L{t, T, = (r - Mt; 0) ■ ■ ■ - fm{t; 0), 

(1.5) M \<f,{t;0-Mt;0\>o. 

Is I — 1 5tGM. 

Let us point out the main difficulties when trying to establish dispersive estimates 
(i.e. the time decay estimates for the U'-L'^ norms) for equation (11. ip . Contrary to 
the energy methods, for dispersive estimates we need to have a good idea about the 
propagators for the Cauchy problem fll.ll) - (ll.2p . Thus, we need to make advances in 
the following two problems: 

• to derive representation formulae for propagators for the Cauchy problem 
(ll.ll) - (ll.2p with time- dependent coefficients. Ideally these propagators would 
be in the form of oscillatory integrals; 

• to analyse the obtained representation formulae for propagators taking into 
account the geometric properties of characteristics which we know should be 
responsible for the time decay rates of U'-L'^ norms of propagators. 

Thus, the aim of the paper is twofold. First, we will present representation formulae 
for propagators for such equations. For this purpose we will develop the asymptotic 
integration method which is a parameter dependent version of the asymptotic inte- 
gration of ordinary differential equations (see e.g. [7|). We will trace the dependence 
on the parameter (which is the frequency in this case) which is essential for further 
investigation. This method, however, will present somewhat surprising results. For 
example, the amplitudes of propagators expressed in this form will have symbolic be- 
haviour of type (0,0) rather than the usual (1,0). Nevertheless, this will be enough 
to carry out the second aim of this part which is the further investigation of the time 
decay properties of the propagators. The price that we will have to pay is that we 
may have to assume additional regularity of the Cauchy data for high frequencies. 
However, this is not so bad because estimates for bounded times already will require 
similar regularity assumptions. 

We will analyse the obtained representations to derive time asymptotics of L^-L^ 
norms of the necessary oscillatory integrals. There are several important differences 
with the case of the wave equation, where level sets of characteristics are nothing else 
but spheres, so one can simply apply the stationary phase method to the obtained os- 
cillatory integrals. Now the critical points may be degenerate so the stationary phase 
method (especially in the parameter depending setting that we have here) does not 
work. In fact, the non- degeneracy of critical points is a rather strong assumption for 
higher order equations, where degeneracy of higher order may easily happen (exam- 
ples of this are e.g. in [2S])- That is why we will allow them to be degenerate of a 
finite order, and the time decay rates will depend on this order. On the other hand, 
van der Corput type estimates that are normally used in place of the stationary phase 
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in such problems are essentially one-dimensional, and they do not take into account 
the geometric properties of characteristics (phases and characteristics do come from 
a hyperbolic equation after all). So, we need to apply a parameter dependent version 
(now time is the parameter) of van der Corput's lemma uniformly in n — 1 directions 
of non-vanishing higher order curvatures. Moreover, this has to be done uniformly 
with respect to the time dependence of the propagators. In fact, we will relate the 
time-decay rates to the Sugimoto's indices of levels sets of characteristic roots of the 
limiting equation, thus establishing a more or less complete picture of perturbation 
properties of dispersive estimates for strictly hyperbolic equations with homogeneous 
symbols. For example, in the case of convex level sets one introduces the convex 
index 7 which is the largest order of tangency of tangent lines to the level sets of 
characteristics of the limiting equation and it turns out that the U'-L'^ norm of the 
corresponding propagator decays as t t . In the case of the second order equations 
one has 7 = 2 and so one recovers that standard rate of decay of the wave equation 
(see in El [TOl [23] ) , and many other known results for the time independent wave type 
second order equations. We also note that such index 7 does not play any role for 
U'-LP estimates, where singularities of the projection from the canonical relations to 
the base space start playing a role (see e.g. survey paper p^). The inclusion of mixed 
derivatives in the symbols may influence the value of 7. Moreover, mixed terms may 
make the analysis more complicated. Already for the second order equations this 
was demonstrated by Hirosawa and Reissig in P], for the problem of the influence of 
oscillations in coefficients. 

In addition, methods introduced in this paper may be applied to the study of 
strictly hyperbolic systems. For example, let A{t, D^) be the first order m x m 
pseudo-differential system, with entries aij{t,^) being homogeneous with respect to 
^ of order one and such that dtaij{-,^) G L^(M) for all C, G K". We consider the 
evolution equation 

(1.6) dtU = iA{t, D^)U, U{0, x) = /(x), x e W. 

Let us assume that system fll.6p is uniformly strictly hyperbolic (see Mizohata [U]), 
i.e. that its characteristics ipkit; k = 1, ■ ■ ■ ,m, are real, and satisfy condition 01.51) . 
Then the strict hyperbolicity implies that we can diagonalise it similar to Lemma [2?T1 
Thus, system (11. 6p splits into m scalar first order equations of the form 

dtVk=ifkit;D^)vk, k = l,...,m, 

for function Vk related to the original vector function U. The condition on the inte- 
grability of time-derivatives of A implies that there is a limiting system A^{t,C,) = 
limf^-too ^(^, with characteristics (^^(0 = linif^±oo <^fc(^; which exist since we 
assume that dtipk{-',0 ^ -^^^(1^)- Therefore, solutions Vk{t,x) can be analysed using 
estimates for oscillatory integrals that we establish in §4. Details of this analysis are 
different from those for the scalar equation (11. II) . especially in the way of keeping 
track of the representation form for the time derivatives of the solution, so we omit 
the analysis of systems from this paper and it will appear elsewhere, together with 
its specific applications and with refinements of the analysis of high frequencies. We 
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will also not discuss the case of oscillations in this paper, but we refer to, for exam- 
ple, the survey [T7], for the overview of the case of the wave equations. The case of 
oscillations in higher order equations will appear elsewhere. 

Thus, in §2 we will discuss the asymptotic integrations of the ordinary differential 
equations corresponding to our problem. Using these implicit representations, we will 
succeed to obtain the asymptotic integrations of fll.ip - (ll.2p . The precise statement 
will be given in §3. 

Let us now give an informal overview of this method. Writing equation (11. ip as a 
system for 

and taking the Fourier transform with respect to x, we can reduce it to the first order 
Cauchy problem 

(1.7) DtU = A{t;OU, U{0) = Uo. 

If we denote ^ 

= Vjis; ds, j = 1,..., m, 
Jo 

a natural candidate for the fundamental matrix for (ll.7p is 

$(t;0 = diag (e^'^i (*'«),.. ■ ,e''^™(*'«)) . 
So, we look for the solution of (I1.7P in the form 

U = ^t;OV, with V{t-0=a{0+£it;0, 

where we want e{t; ^) to decay as t ^ ±oo. It can be checked that there is a matrix 
y4o(t; such that = Aq^ and we get 

DtU = Dt{(!>V) = {Dt^)V + ^DtV = AqU + ^DtE. 

Thus, U becomes the solution of (11.71) if we choose V and e such that 

(1.8) DtV = Dte = ^-\A - Aq)^V. 

In §2 we will show that, in fact, there exists a global- in-time solution V of equation 
(II. 8p of the required form V = a+£. Moreover, £ satisfies the property that £{t;^) 
as t ^ ±oo for all ^ 7^ 0. In addition, we will show the decay orders of both a{-) and 
e{t] ■) and their derivatives. This will lead to an oscillatory integral representation of 
solution u{t, x) of (II. ip of the form 

m— 1 m 

(1.9) u{t,x) = J2Y.^^' 

k=0 j=l 

with amplitudes _^{^) , el _^(t; C.) of the form of a and £ above. In fact. Theorem 
13.11 will also yield a similar representation for the derivatives of u{t,x) with respect 
to time. The main difference with equations with time independent coefficients here 
is that the amplitudes al^{^) and el^(t;^) will have the symbohc behavior of the 
type (0,0) rather than the type (1,0) usual for equations with constant coefficients. 
Indeed, such choice of phases globally as 'dj(t;^) introduces low order errors in the 
equation if we formally substitute (11.91) into (ll.ip and as we know the lower order 
terms may change the time decay properties in an essential way (this is especially 
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apparent for Schrodinger equations, but is also true in the hyperbolic case). Thus, 
the error should be somehow accounted for and the behaviour of amplitudes takes 
care of this. In any case, since we know that the needed regularity of data comes from 
other parts of the time-frequency phase space, we are still able to get the same time 
decay rate under an additional regularity assumption in the high frequency zone. So 
this difference does not change the final result in a big way. 

Thus, in the second part of the paper we will use representation (11.91) to derive the 
time decay of the L^-L'^ norms of m, which in turn leads to Strichartz estimates and 
to well-posedness results for the corresponding semilinear equations in a rather (by 
now) standard way, so we will derive the dispersive estimates and will omit the details 
of the further standard analysis. In addition, in §4 we will present estimates for more 
general oscillatory integrals. Such estimates may be used not only in the application 
to the problem we are considering in this paper but in a wider range of applications. 
The estimates will rely on estimates for parameter dependent oscillatory integrals 
developed in [20]. We state such result here in Theorem 14.51 However, the meaning 
of the parameter is different in our setting. Thus, in our problem here time acts as a 
parameter while in problems for hyperbolic equations with constant coefficients but 
with lower order terms considered in [22] the phase functions were not homogeneous 
and their non-homogeneous contributions were considered to be a parameter from 
the point of view of the perturbation theory. In principle, it should be possible to 
combine problems with time-dependent coefficients with those with lower order terms 
but this will be a subject of another paper - here we have an advantage of making 
more use of the homogeneity of the symbols and hence also of phases, considerably 
simplifying some arguments. The obtained results can be applied to the global in 
time well-posedness problems of Kirchhoff equations of high orders and of Kirchhoff 
systems. Such applications will be addressed elsewhere. 

Let V'fc (0 = liHit^itoo '^k{t'i i) be the limits of characteristic roots as will be shown 
to exist in Proposition 12.31 Let us introduce the convex and non-convex Sugimoto 
indices for the level sets of these functions. In the time independent setting these 
indices have been introduced by Sugimoto in [2U These indices will determine 
the decay rate of propagators for large frequencies. 

Let G C°°(M'^\0) be a homogeneous of order one function and let S,^ = G M" : 
ip{S,) = 1} be its level set. Suppose first that T,^ is convex. We define the convex 
Sugimoto index 7(S^) of by 

(1.10) l{^<fi) '■= sup sup7(S<^;cr, P) , 

o-es^ p 

where P is a plane containing the normal to Tj^ at a and 7(S(^; a, P) denotes the 
order of the contact between the line f} P (where is the tangent plane at a), 
and the curve T,^ fl P. 

In the case when the level set is not convex, we get a weaker result based on the 
van der Corput lemma. In this case we use the non-convex Sugimoto index 'JqCE^) of 
Tj^ which we define as 

(1.11) 7o(S<^) := sup inf 7(S^; a, P) , 
where P and a are the same as in the convex case. 
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We note that for the second order equations we have 7 = 70 = 2 and the following 
theorem covers the case of the wave equation as a special case, also improving the 
corresponding result in [T2]. We use the notation L^{W^) for the standard Sobolev 
space with s derivatives over L^(M"), and by L^(]R") we denote its homogeneous 
version. The result on the dispersive estimates that we will prove among other things, 
is as follows: 

Theorem 1.1. Assume fll.3p -( fT31) . Then the solution u(t,x) of (11. ip satisfies the 
following estimates: 

(i) Suppose that the set 

= G : if^iO = 1} 
is convex for all k = 1, . . . ,m, and set 7 = max 7(S ±). In addition, suppose that 

k=l,...,m 

(1 + \t\ya'^^j G L^(R) for 1 < r < [{n - l)/7] + 1, and for all u,j with |z/| + j = m. 
Let 1 < p <2 < q < +00 and ^ + ^ = 1. Then for alltEM. we have the estimate 



1.12) pjD>(t,-)l|L.(R")< 



m— 1 



c{i + \t\)-^ihl) (llMk.^^^,^,^,_,(M«) + IIMkr,„.„(M^)) 

k=0 



where Np = (n — 
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/ = 0, . . . , m — 1, and a is any multi- 
index. 

(ii) Suppose that S^^i is non-convex for some k = 1, . . . ,m, and let us set 70 = 
max 7o(S ±). In addition, suppose that {l + \t\)al , G L^{Wj for all u, j with\v\+j = 

k=l,...,m " 

m. Let l<p<2<q< +00 and ^ + ^ = 1. Then for allt we have the estimate 

P^D>(t,-)|U.(M")< 

m— 1 

C(l + \t\)-i(H) J2 (llAlli.„.,„.,_.,r, + IIMIif,.,-.<«")) • 

k=0 

where Np = (^n — ^ -\- 1^ (p ~ g) ' ^ = 0) • • • 5 '"^ ~ 1; ^'^^ is any multi-index. 

Remark 1.2. The way we formulate the estimates in Theorem \l.l\ is to unify different 
estimates for different parts of the solution. This may explain the appearance of two 
norms in the right hand side of (11.120 . for example, to account for both small and 
large frequencies. The much more precise estimates are possible and they are stated 
in Theorem \4.1C\ 

Let us now make only a few short remarks to compare our results with what is 
known for m = 2. In the constant coefficient case and p = 2 the estimates coincide 
with those for constant coefficient equations considered in [21]. Also, in the case of 
constant coefficients, the Sobolev index Np in the estimate (I1.12p can be improved 
since there is no addition of the integer part in its definition in this case. For the 
detailed overview of constant coefficients case we refer to |22], and results in this 
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direction for non-constant coefficients were announced in [H] and the detailed proofs 
will appear elsewhere. 

Further, in the time-dependent case, it can be already noted that the statement 
of Theorem 11.11 goes beyond results available in certain energy classes. For example, 
in the often considered case of the time-dependent wave equation (so that m = 2, 
e.g. [8], [13], [16], [17], etc.) one obtains the estimate for \\Dtu(t, ■)\\li{R''i) and 
||V'u(t, •)||l9(r") only, and not for solution u itself. Moreover, the use of homogeneous 
Sobolev spaces in fll.l2p allows to gain more information in the low frequency region. 
At the same time, also already for the case m = 2, we make the assumption on only 
one derivative of the coefficients a^j, which is another improvement compared with 
the known literature. This improvement will be crucial in dealing with applications 
to Kirchhoff equations. 

We will denote {x) = ^1 + Constants may change from formula to formula, 
although they are usually denoted by the same letter. 

The authors thank Jens Wirth for remarks leading to an improvement of the man- 
uscript. 

2. Asymptotic integration of ODE 

In this section we will construct the asymptotic integration of the ordinary differ- 
ential equation. By applying the Fourier transform on to (11. ip . we get 

m 

(2.1) D^v + Y,h,{t-ODT-'v = 0, 

where 

(note that there is a slight change of the meaning of j here compared to (II. ip ). This 
is the ordinary differential equation of homogeneous m*^ order with the parameter 
^ = (,^1, . . . , ^n)- As usual, the strict hyperbolicity (ll.4l) - (ll.5p means that the charac- 
teristic roots of (12.11) are real and distinct. We denote them hj (pi(t; ^) , . . . , (pmit; C,). 
Notice that each ipiit;^) has a homogeneous degree of order one with respect to ^. 
In this section we will develop an asymptotic integration of the equation (12. ip . Let 
us start by writing (12. ip as the first order system. In (12. ip we put for brevity 

and denote 

v, = \^r-^-'Dlv, j=0,l,---,m-l. 

It is easy to see that 

DtVj = \^\vj+i, j = 0, 1, • • • , m - 2, 
holds. Then (12. ip can be written as 

m— 1 
j=0 
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Hence, if we put 



n{t,o 



/ 1 ... \ 

■•. 

; ■•. ■•. 1 

then (12.11) can be written by 

(2.2) D,v = n{t,Cmv, 

where v = '^{vq, f i, . . . , Wm-i)- 
We will use the following lemma. 

Lemma 2.1 ([15j Proposition 6.4). Assume (11.31) -( |T31) . Then there exists a matrix 
Af = A/'(t; ^) of homogeneous degree satisfying the following properties: 

(i) Afnit,C) = VAf, where 

V = v{t- = diag {^i(t; e/iei), • • • , ^m{t; i/m , 

(ii) inf IdetMt;^! > 0, 

geiR"\o,teK 

(iii) A/'(t;0 ^n ^ ^ 0, CMn t G M and dtAf{t;^) belongs to L^{R) for each 

We will first derive the energy estimates. 

Lemma 2.2. Assume (11.31) -f fT75|) . Let v = v{t]C,) be a general solution of (12.21) . 
Then, for all t G M, we have 

(2.3) \v{t;Of < C|t)(0;O|V-- 

Proof. Multiplying (12.21) by A/" = A/'(t; .^) from Lemma 12.11 we get 

Dt{Afv) - .Un\^\v - {Dt^)v = 0. 
Putting Mv = w, we have 

(2.4) Dtw - V\^\w - {DtU)v = 0, 
since AfTi = T>M by Lemma 12.11 This implies that 

dt\w\^ = 2Re {dtW ■ w) = 2Re {iV\^\w ■ w) + 2Re {i{DtM)v ■ w) . 



Taking account that iT>\(^\ = —iT)\^\ and T> is real and diagonal, we have 

Re{iV\^\w-w) = 0, 

hence, 

(2.5) dt\w\'^ <2\\dtAf\\\v\\w\. 

Here, \v\ and \w\ are equivalent to each other. Indeed, there exists Ci,C2 > such 
that Ci\v\ < \w\ < C2\v\ on account of Lemma 12. 1[ Thus, integrating (12. 5p . we arrive 
at 

\vit; P < C (^viO; OP + \\dtAf{s; II Hs; P ds^ . 

bmce dtAf G L\R) by Lemma O (iii), we conclude from Gronwall's lemma that 
(12. 3p is true. The proof of Lemma [2.21 is complete. □ 
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As a consequence of (12 ■4p in the proof of Lemma 12.21 we have derived 

(2.6) Dtw - V\^\w - {DtU)U-^w = 0, 

where we put w = Afv. We can expect that the solution of (12. 6p is asymptotic to the 
solution of 

(2.7) Ay = my- 

Let $(t;^) be the fundamental matrix of (12. 7p . i.e., 

$(t;0 = diag{e''^^(*'«),--- , e''^™(*'«)} , 

where we put 

^i(^;0=/ <^j{s;Ods, j = l,...,m. 
Jo 

Let us first analyse certain basic properties of characteristic roots ipk{t', of (ll.4p . 

Proposition 2.3. Let the operator L{t, Df, D^) satisfy the properties (II. 4p - (II. 50 . 

Then each dtfkit', 0? ^ — ^, ■ ■ ■ ,tTi, is homogeneous of order one in ^, and there exist 
a constant C > such that 

(2.8) \dtMt;0\<C\^\ for all ^eW, t G R, k = l,...,m. 

Moreover, if alj^-) G L^(R) for all h',j, then we have also dtfk{-',0 ^ -^^"'^(R) for all 
^ G M". Furthermore, there exist functions (f^ G C°°(]R"\0), homogeneous of order 
one, such that 

(2.9) d^Mt; ^ d^^kiO o^s t ^ ±oo, 

for all ^ G M*^, all a, and k = 1, . . . ,m. Finally, we have the following formula for 
the derivatives of characteristic roots: 

(2.10) dtMt;0 = - E <,(i)v^fc(^;0'rn(^^(^'^)-^-(^;0)"'- 

\u\+j=m r=ik 

Proof. Let us show first that (fk{t', is bounded with respect to t G M, i.e., 

(2.11) \^kit;^)\<C\^\, for all ^ G M", t G R, A; = l,...,m. 

We will use the fact that ipk{t', roots of the polynomial L of the form 
L(t, r, = + ci (t, Or"^-' + ■ ■ ■ + c„(t, 

with \cj(t,^)\ < M|^p, for some M > 1. Suppose that one of its roots r satisfies 
|r(t,OI >2M|e|. Then 



\ci{t,^\) 


|Cm(t,OI 




r 




r 


m 



> ItI'"! 1 

>2A/|{r(l-i-^-...-^4^)>0, 

hence |r(t,0| < 2M|^| for all ^ G M". Thus we establish fl2llil . 
Differentiating (II. 4p with respect to t, we get 

^^^^|^= E <.(^K'r = -E5t^.(t;on(^-^^(^;0)- 

\i>\-\-j=m k=l r^k 
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Setting T = ipkit; we obtain 

(2.12) dtMt;Oll^Mt;0-Mt;0) = - E <At)Mt;0'C, 

Vy^k \u\+j=in 

implying fl2.10p . Now, using (11. 5p . (12.111) . and the assumption that a^j(-) G L^(]R) 
for all z/,j, we conclude that ( j2^ holds and dtipk{-]0 ^ L^{^) for all ^ G M" and 
k = 1, . . . ,m. The homogeneity of order one of dt<fk{t', ^) is an immediate consequence 
of (12.121) and its derivatives. 
Finally, setting 

/•itoo 

Jo 

we get (12. 9p with a = 0. Differentiating this equality with respect to ^, we get (12. 9p 
for all a. The proof is complete. □ 

We note that under the assumptions of Proposition 12.31 the coefficients a^j {t) of 
the operator L(t, Dt, D^) in (II. ip have limits a^j as t — > ±oo, namely 

/•itoo 

<i = «..i(0)+ / </t)dt. 
Jo 

Functions (0 characteristics of the limiting strictly hyperbolic operator 

(2.13) L^iD,,D,)u^Dru+ Yl <^D:Diu, 

\i/\-\-j=7n 

J < m ~ 1 

and their geometric properties are responsible for the time decay of solutions to 
the Cauchy problems for both operators L{t, Dt, D^) and L^{Di,, D^). This will be 
analysed in §4. We also note that since operator L^{Dt, D^) has constant coefficients 
its solution can be represented as a sum of oscillatory integrals in the standard way. 
The dependence of coefficients of L{t, Df, D^.) on time brings corrections to the phases 
and amplitudes of this representation. 

Next we make the representation formulae of solutions for our equation. The fol- 
lowing proposition is known as Levinson's lemma (see Coddington and Levinson [6]) 
in the theory of ordinary differential equations, the new feature here is the additional 
dependence on ^. For the convenience of the readers, we shall prove it along the 
method of Ascoli [2] and Wintner [27] (cf. Hartman \7]). 

Proposition 2.4. Assume (ll.3p - (ll.5p . Then, for every nontrivial solution v{t;^) 
of (12. ip . there exist vectors of -amplitude functions C(±{C,) and error functions 
e±(t; ^) such that 

(2.14) vit; = ^fit■, 0-'^t; (a±(0 + £±{t; 0) , 0, 
where 

(2.15) £±{t',0 ~* f'^'i" ^^2/ fixed ^ 7^ as t ±oo. 
Furthermore we have 



(2.16) 



Ae±(t;0 = C(t;0 (a±(0 +e±(i; 0) , 
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where C(t; ^) belongs to L^{M.) in t, and has the following form: 

(2.17) C{t- = 0'\D,Afit; OWit; O^'Ht; 0- 

Proof. We can expect that every solution w = w{t; C,) of (12.61) is asymptotic to some 
solution y = y{t] ^) of (12. 7p . If we perform the Wronskian transform z = <l>(t; ^)~^w, 
then the system (12. 6p reduces to a system DfZ = C{t]^)z, where C{t;^) is given by 

(imD- 

We will now prove that (I2.14I) - (I2.15I) hold for every nontrivial solution v = v{t] ^). 
It follows from Lemma 12.21 that 

k(t;OI < mt;0-'\\\w{t;0\<c\w{t;0\<ci\w{0;0\ 

for all t G R and some constant Ci. Using this bound and equation DfZ = C(t;^)z, 
we have 

\Dtz{t;0\ < \\C{t-m \z{t;0\ < c^\w{0■mC{t■Ol 

hence Dtz{-;^) G L^(R) on account of C{-;C) G L^(]R). Thus {z(t;^)}teR is a conver- 
gent function, and there exists 

lim 2(t;0 =:a±(0- 

If we set 

(2.18) e±{t-0 = z{U)-c^±{0, 

then z{t]^) can be written as z(t;^) = a±{^) +£±(t;^) for t ^ 0, and further, 
£±(^;0 decays as t — >• ±oo for any fixed ^ 7^ 0, which proves (12.151) . Since w(t;^) = 
$(t; ^)z(t; ^), we get the formula (I2.14p . Finally, differentiating (12.180 with respect 
to t and using the equation DfZ = C(t; ^)z, we get (I2.16p . The proof of Proposition 
12.41 is now complete. □ 

Finally, we will need the estimates for higher order derivatives of C{t]^) appearing 
in Proposition | 



Lemma 2.5. Assume (11.30 - (11.50 . Then the fj}^ derivatives for each entry Cjkit;^,) 
ofC{t]0 satisfy 

(2.19) |9f c,fc(t; 01 < c\i\-\^\{l + \t\t\^{t), j,k = l,...,m, 
for \fi\ > 1 and \C,\ > 1, where 

(2.20) m= E KM- 

\u\+j=m 

J < m — 1 

For < 10 < 1 and > 1, we have 

(2.21) \d^c,k{t; 0| < c|0-l^l(l + Itl)l'^l^(t), j,k = l,...,m. 

Moreover, assume that (1 + \t\y^'\alJ{t) G L\m) for all and for some (X with 
> 1. Then we have d^Cjk{t;C,) G -/^-"^(M). 
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Proof. Since fj{t] ^) is homogeneous of order one, we have 

\d^¥>At-^ 01 < c|er'^l+' for e e \ 0, J = 1, . . . , m, 

and hence, 

r\t\ 

\d^Mt;0\< / \d^vAs;0\ds<c\tm-^^^+\ foreGM"\0, j = l,...,m. 
Thus we get, for every multi-index /x. 

Now let us go back to fl2.17l) . It follows from (12.101) that dtAf is represented by 
alj. Hence, using (12.221) and differentiating (12.171) with respect to ^, we conclude 
that jj}^ derivative of Cjk{t; with respect to ^ is bounded by |0~'^'(1 + |t|)'^'\E'(t) for 
10 > 1, where "^(t) is given in (12.201) . d^Cjk{t] G L^(M) follows from the assumption 
(1 + \t\y^^alj{t) G L^(M) for all uj. The proof of Lemma [221 is complete. □ 

3. Representation of solution 

In this section we will establish the representation formulae for solutions of the 
Cauchy problem (12.11) in the form of the oscillatory integrals. Let 

V{t;0='^{vo{t-0,.-.,Vm^i{t;0) 

be the fundamental matrix of (12. 6p . This means that V^(0;O = I- Then it follows 
from Proposition 12.41 that each Vj{t;^) can be represented by 

(3.1) v,{t; = O^'Ht; K,±(0 + e.,±(i; 0) • 

Let u{t,x) be the solution to (11.11) with the Cauchy data D^u{0,x) = fk{x). Put 

uit, = ii^r-'uit, 0, • • • , \^r-'-'Diu{t, 0, . . . , DrMt, of ■ 

Then we can write u(t,0 = V(t;0'^{^yOi ^^us we arrive at 

m m—l 

(3.2) \^r-'-'Dlu{t,0 = 5^ 5^e^'^(*^«)n'^(t;0 K±(0 + 4±('^; 0) lO'^-'-'M), 

j=l k=0 

where n''^{t] is the entry of Afit; 0^^- 



.,771—1 
j = l,...,m 



Summarizing the above argument, we have the representation formulae of (II. ip - 

Theorem 3.1. Assume (11.31) - (11.51) . Then there exists al. and e-'f, _^{t] ^) such that 
the solution u{t,x) of our problem (II. II) - (11.21) is represented by 



m m—l 



D\u{t,x) = EE-^"' p(*^«)n'^(t;0 K±(0 \0'-%{0 

j=l k=0 
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for I = 0, . . . ,m — 1, where 

r+oo 

J\t\ 

and ^(t) is defined by fl2.20p from Lemma \KR For the higher order derivatives of 
amplitude functions, we have, for > 1, 

< c, \dleUt;i)\ < ce/<l'(l+^)""*(^)'^^ |^| > 1, 

// in addition to (11.31) -( fT3i) . we further assume that (1 + |t|)l'^l<^j-(t) G Li(M) for 
some fi with > 1, and for all z/, j, then the bound for each d^e^j_{t;^) is uniform 
in t. 

Proof of Theorem \3.1[ We must determine the precise growth order of and 
with respect to ^. 

Lemma 3.2. Assume (11.31) - fll.Sp . Then there exists a constant c > such that, for 
j = 1, . . . ,m and k = 0, . . . ,m — 1, 

/■+0O 

\al^{0\<c, |4±(t;0| <c / ^{s)ds. 

J\t\ 

Proof. Let us go back to (13.11) . Then we have, by using Lemma [2.2[ 

(3.3) K-,±(0+£,,±(t;OI < \^t-iY^M{t-i)v,{t,o\<co 

for j = 1, . . . , m. Since £j,±{t] ^) decays to 0, it follows that 

\Oij,±{0\ < Co. 

On the other hand, (12.161) is equivalent to the following integral equation: 

r+oo 

£j,±{t; = M C{s; {£j,±{s; + otj,±{0) ds. 
J\t\ 

Thus combining this equation and (13.31) . we get 

f+oo r+oo 

|e,,±(t;OI <Co / \\C{s-mds<c / ^{s)ds. 
J\t\ J\t\ 

The proof of Lemma 13.21 is finished. □ 

We need the estimates of higher order derivatives of amplitude functions. 
Lemma 3.3. Assume (11.31) - (11.51) . Then we have, for > 1, 

|a[4±(o| <cier'^', o<iei<i, 

(3.4) |'9f4,±(t;0| < ce/o'(l+^)""*(^)'^^ \i\ > 1, 

(3.5) |9f4±(t;0| < ce^o"(i+^)""*(^)^^|erl'^l, < |e| < 1. 
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In addition to fll.3p - fll.5p . if we assume that (1 + e L\R) for some n 

with > 1, then (13.40 -f l375l) is uniform in t. 

Proof Putting 

£m— 1,± 

we see that the matrix Q{t,C,) satisfies 

DtQit,0=Cit;OQit,0 

with the initial data 

Q(0, = m O"W(0; K(0; 0, ■ ■ ■ ,«™_i(0; 0) = A/'(0; 0- 

Then it follows from the theory of ordinary differential equations that Q{t,C,) can be 
written by Picard series: 

(3.6) Q{t,0 = 

(i + ^l Cin; dn + ^ Cin; dn £ C{t2, i) dT2 + -- A/'(0; 0- 
We note from Lemma [2.51 that 

'c(l + |t|)HvI/(t), > 1, 



(3.7) ||afC(t;0|| < 

where 



c(i + |t|)H^(t)ieri^i o<iei<i. 



j-\-\y\=m 

j<m — l 

Differentiating (13.61) with respect to C,, we have, by using (13. 7p . 
(3.8) |af (a,,±(0 < ce/o"(i+^)""*(^)'^^ 

for all t G M, |.^| > 1, > 1 and j = 0, . . . , m — 1, where we have used the following: 
Fact. Let f{t) G C(R). Then 

efsfi-)'^- = 1+ l"f(^r,)dn+ ff{n)dn H f{r2)dT2 + ■ ■ ■ . 

J s J s J s 

Since dtd^Q{t,0 = {dtd^£l.^^{t;0), we combine DtQ = CQ and (13. 8p to deduce 
that dtd^£''^{0;C,) exists and is uniformly bounded in |^| > 1, which ensures the 
existence of 9^e-j_(0;^). Using again (13.81) with t = 0, we conclude 

In a similar way, we get the bound for 9^Q!j,±(0 low frequency part < |,^| < 1. 
If we combine these estimates with (13.81) . we have the bound for d^ej^±{t; ^): 

The proof of Lemma 13.31 is complete. □ 
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Completion of the proof of Theorem \3.1[ The estimates of the amphtude and error 
functions have been derived in Lemmas I3.2H3.3I The proof of Theorem 13.11 is now 
finished. □ 



4. Estimates for oscillatory integrals; Proof of Theorem 11.11 

The aim of this section is to estabhsh time decay estimates for U'-L'^ norms of 
propagators for the Cauchy problem (11.11) . which gives the proof of Theorem II. 1[ 
The analysis of high frequencies will give estimates dependent on the geometry of the 
level sets of characteristic roots of the equation. For small frequencies estimates are 
independent of the geometry of the level set and are given by Proposition 14.11 below. 
We recall that Theorem 13.11 assures in particular that the solution to the Cauchy 
problem (11. ip is of the form 

m— 1 m 

(4.1) u{t,x) = Y,Y.^' 

k=0 j=l 

The following Proposition 14.11 provides the time decay estimate for small frequencies 
for each of the terms in this sum. To simplify the notation, we formulate it in a more 
general form for general oscillatory integrals of the form 

TJ{x)= [ e^(^-«+^(*^«))a(t,0/(0^e. 

In the analysis of oscillatory integrals in the sum (14. ip we will actually make time- 
dependent cut-offs and analyse separately different ranges of frequencies. We can 
obtain the following proposition for small frequencies |^| < t~^. Higher frequencies 
1^1 > t~^ will be analysed later. Thus, we introduce a cut-off function of the form 
^'((l + |t|)0 for some G Co"^(M") such that V^(0 = 1 for |^| < i, and for 
1^1 > 1- We recall that we use the notation L^iW^) for the homogeneous Sobolev 
space W^{W). 

Proposition 4.1. Let T^, t G M, be an operator defined by 

Ttf{x) = [ e^(-«+''(*^«))V^((l + m)a{t, 0/(0 d^, 

where ^{t;C,) is real valued, positively homogeneous of order one in ^. Assume that 
the amplitude a(t, satisfies 

(4.2) l01«(i,OI<C 

for some k G M, and for all t & K. and all ^ G suppip{{l + |t|)0- Let 1 < p < 2 < 
q < +00 be such that ^ + ^ = 1. Then for t eM. we have the estimate 

(4.3) ||T,/|U,(M.) < C(l + JR.), 
where constant C depends on n,p,q and the norm || l^l^o-Hioo . 
Proof. We easily obtain 



K±(0 + 4,±(^;0)e''^(*^^)/.(0 



t ^ 0. 
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by the Plancherel identity. Thus (14. 3 p would follow by analytic interpolation from an 
estimate: 



(4-4) \\Ttf\\L--{R") <C{1 + \t\) WJWllJRny 

In fact, since 011/(01 — C*]] lO'^c^lU"" II/IIl?_^(R")) "we can estimate 

iiTt/iu^(Mn) < / ia(t,oii7(OMe 

<c( [ dA\\f\\L^ju^)<Cil + \t\)-^ 



I5l<(l+|i|)-^ 

for all t G M. This proves (14. 4p . The proof of Proposition 14.11 is complete. □ 

We now turn to the analysis of larger frequencies. The following proposition pro- 
vides the necessary background to obtain the time decay estimate for large frequencies 
for each of the terms in the sum (14. ip . In fact, it will be used for frequencies |0 ^ 1 
but will be formulated here in a slightly more general form. The relation with the 
sum (14. ip and a refinement for low frequencies will be also made in Proposition 14.61 

Proposition 4.2. Let Tt, t > 0, be an operator defined by 

T,f{x)= I e^(^-«+''(*^«))a(t,0/(Orfe, 



where i?(t;0 '"^a/ valued, continuous in t, smooth in C, & M"\0, homogeneous of 
order one in ^. Assume that the set 

= G M"\0 : ^(0 = 1} 

is strictly convex and let 7 = 7(S^) be the convex Sugimoto index of S<^, as defined 
in (ll.lOp . Suppose that 

W{t-i)\<C{l + tm for all t>0, eeK", 

and that there is some ip G C°°(]R"\0), > 0, such that 

(4.5) t-^d'^^it] ^ d'^^{i) as t^ 00, for all ^ G R"\0, |a| < 7. 
Assume also that the amplitude a(t,0 satisfies 

(4.6) \d^a{t,0\<Ca for all |a| < [(n - l)/7] + 1. 

Let l<p<2<q<oo be such that ^ + ^ = 1. Then for t > we have the estimate 



P Q 

n-1 f 1 1 ' 



(4.7) mfh^^n) < 



where N„ = [n — + 



p \ 7 



n-1 



+ 1 



1 _ 1 
p q 



The number of derivatives Np = — + + 1 j — ^ j required for the 

estimate (14. 7p is determined by the fact that the amplitude a(t, i^i (14. 6p is in the 
symbol class Sqq rather than the usual S^q. In fact, if a(t, satisfies inequalities 

I 11 — 1 

(4.8) |5f a(i,OI < C'c.(0" for all lal < +1, 

« L 7 
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then we can take the Sobolev index Np = {n — — for the estimate (14. 7p to 

hold. However, the method of asymptotic integration and the statement of Theorem 
13.11 forces us to assume (14.61) rather than (14.81) . 

Let us now discuss other assumptions we make in this proposition from the point 
of view of the original Cauchy problem (II. ip . We recall from (12.130 that functions 
V^fc (0 characteristics of the limiting strictly hyperbolic operator 

(4.9) L^iD,,D,)u^Dru+ aJ^-D^D^w, 

\i/\-\-j=7n 

< m — 1 

and their geometric properties are responsible for the time decay of solutions to 
the Cauchy problems for both operators L{t, Dt, D^) and L^{Dt, D^)- The fact that 
(p^{^) are characteristics of (14. 9p . implies that they are real analytic for ,^ 7^ and that 
we have the following statement, which was established for operators with constant 
coefficients by Sugimoto 

Proposition 4.3. Letipk{C), k = 1, . . . ,m, be characteristics of operator (14.91) . or- 
dered by ipi{C,) > (p2{0 > ■ ■ ■ > V^m(0 fo'i^ 7^ 0. Suppose that all the Hessians (pliO 
are semi- definite for ^ 7^ 0. Then there exists a polynomial a;(^) of order one such 
that ^pm/2{0 > «(0 > ^m/2+1 {iffn is even)j)r a{i) = '^(m+i)/2iO {^f ^ odd). 
Moreover, the hypersurfaces T.^ = {C, E M"; ip^. = ±1} with (pk{0 = V'fclO ~ 
(fc 7^ (m + l)/2) are convex and 7(Sfc) < 2[m/2]. 

In particular, in our arguments we can replace ipk by (fk since the addition of 
a linear function does not change the decay rate nor the index 7(S(pj,). This also 
ensures that the limiting phase in Proposition 14.21 may be taken to be strictly 
positive. Indeed, it can be taken to be nonzero, and if it is strictly negative we 
simply replace ip by —ip. Moreover, the assumption that is strictly convex in 
Proposition 14.21 can be replaced by the assumption that it is only convex. Indeed, 
since ip{^) is a characteristic root of (14.91) . it is real analytic for ,^ 7^ 0. Then, the 
convexity, the real analyticity and the compactness imply that it is actually strictly 
convex. In particular, it also implies that 7(S<^) is finite and even. 

In the case when the level set T,^ in Proposition 14.21 is not convex, we get a weaker 
result based on the one- dimensional van der Corput lemma. In this case we use the 
non-convex Sugimoto index of which was defined in (II. lip in the introduction. 

Proposition 4.4. Let Tt, t > 0, be an operator defined by 

Ttf{x)= [ e^(-«+'^(*^«))a(t,0/(Orfe, 

where ^{t;C,) is real valued, continuous in t, smooth in C, E R"\0, homogeneous of 
order one in ^. Let 70 = 7o(S(^) be the non-convex Sugimoto index of the level 
surface = G ]R"\0 : ip{^) = 1}. Suppose that 

mt-0\<C{l + tm for all t>0, eeR^ 

and that there is some ip G C°°(M"\0), ip > 0, such that 

t-^d^^{t; d^ifiiO as t^oo, for all ^ G M", |a| < 70. 
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Assume also that the amplitude a(t,^) satisfies 

(4.10) \d^a{t,^)\ <Ca for all \a\ < 1. 

Let l<p<2<q< +oo be such that - + - = 1. Then for t > we have the estimate 



p <? 



(4.11) \\Ttf\\L^iUn^<Ct ^^||/|U5^^(Mn), 

whereN,= (n - + 1^ (p " ' 

We will first prove Proposition 14.21 and then indicate the changes necessary for the 
proof of Proposition | 



Proof of Proposition^?^ First we observe that ||Tt/||/^2 < C||/||j;^2 by the Plancherel 
identity. In order to simplify the proof somewhat, we will absorb the Sobolev index 
Np into the amplitude a(t,^), so that we will estimate ||Tt/||i<j in terms of 
and not in terms of 1 1/| . Thus, instead of (14.61) . from now on we will assume that 

the amplitude a{t,S,) satisfies 

\d1a{t,i)\<C^{i)-' 

for all \a\ < [(n — l)/7] + 1 and k = Np, so that estimate (14. 7p would follow by 
interpolation from the estimate 

(4.12) ||Ti/|Uoo(Kn) < Ct-^ 



and where we take k = N = n — + 



n-l 



1. Note that since we assume that 



7 

the amplitude is bounded for small frequencies, we can work with standard Sobolev 
spaces here. By using Besov spaces, we can microlocahse the desired estimate to 
discs in the frequency space. Indeed, let {$j}°2,Q be the Littlewood-Paley partition 
of unity, and let 

(oo 
j=Q 

be the norm of the Besov space -Bpg- Then, because of the continuous embeddings 
LP C for 1 < p < 2, and B^^ C for 2 < g < +00 (see 0), it is sufficient 
to prove the uniform estimate for the operators with amplitudes a(t, ^)$j(^). Let us 
denote 

dit,o=t-'^it,o, 

so that by the assumption we have '&{t;C,) —>■ y^iC,) as t ^ 00. Now, writing 




$,(0 = $,(0^l 



2i 



with some function G C^(0, 00), we may prove the uniform estimate for operators 
with amplitudes a(t, ^)\I^ ^^^j^ j . Such choice of \I' is possible due to our assumption 
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that ^{t; ^) (p{^) as t ^ oo, and we restrict the analysis for large enough t. Let 

(4.13) I{t, X) = £ e^^-^^'^^-'^Mt, 0^ 

be the kernel of the corresponding operator. Since we easily have the L^-L^ estimate 
by the Plancherel identity, by analytic interpolation we only need to prove the L^- 
L°° case of fl4.12p . In turn, this follows from the estimate |/(t,a;)| < Ct ~, with 
constant C independent of j. 

Let K, G C^(R"') be supported in a ball with some radius r > centred at the 
origin. We split the integral in 

I{t,x) = Ii{t,x) + h{t,x) 

+ j^^ e^(-^+'^(*^^))a(t, 0(1 - K) (t- V + t- V,^(t; 0) ^ j d^. 

n — 1 I / \ I ; 

We can easily see that |/2(t, x)| < Ct ~ . In fact, we can show |/2(t,a;)| < Ct~ for 
I = [{n — l)/7] + 1 and then the required estimate simply follows since I > {n — l)/7. 
Indeed, on the support of 1 — k, we have \x + V ^'d{t] ^)\ > rt > 0. Thus, integrating 

by parts with operator P = ^|^^v^^^{*.g}|2 ■ V^, we get 

(4.14) l2{t,x) = 



a(t,0(l-«:) {t-'x + t-'V^^it-O)"^ 



2i 



di. 



Using the fact that \d'^^{t;C)\ < C(l + 1)|0^"'"', we readily observe from fOij) that 
the required estimate |/2(t,x)| < Ct~^ holds. Here we also used the condition (14. 6 p 
which assures that we can perform the integration by parts [{n — l)/7] + 1 times. 

Now we will turn to estimating Ii{t,x). Recall that ^{t;C,) = t^^i}{t;C,) and 
■^(t; ~^ viO as t ^ oo. Let us denote 

S* = {eeM":^9(t;0 = l}- 

It can be readily checked that 7(S*) 7(S<^) = 7 as t — oo. So we can restrict our 
attention to t large enough for which we have 7(S*) = 7. By rotation, we can always 
microlocalise in some narrow cone around e„ = (0, . . . , 0, 1) and in this cone we can 
parameterise 

^' = {iy,ht{y)):yeU} 

for some open U C M"^^. In other words, we have d{t\y,ht{y)) = 1, and it follows 
that ht is smooth and Vht : U Vht{U) C M*^"^ is a homeomorphism. The function 
ht is concave if S* is convex. We claim that 

(4.15) \d^ht{y)\ < Ca, for all y G U and large enough t. 
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Indeed, let us look at |a| = 1 first. From {}{t] y, ht{y)) = 1 we get that 

From homogeneity we have \W^'d\ < C, so also {Vyi^l < C. By Euler's identity we 
have 



(4.16) 



d^^d{t] e„) = §{t] e„) v^(e„) > as t ^ oo, 



so we have \d^^'d\ > c > since we are in a narrow cone around e„. From this it 
follows that \Wyht{y)\ < C for all |/ G t/ and t large enough. A similar argument 
proves the boundedness of higher order derivatives in (I4.15P . 

Now, let us turn to analyse the structure of the sets S*. We have the Gauss map 



Vc^(t;C) 



5n— 1 



|Vc^(t;C)l 

and for x = {x', x„) G R"^^ x M near the point — V('i9(t; e„) we define G f/ by 

{zt, ht{zt)) = u'^{-x/\x\). 
Then {—'Vyht{y), 1) is normal to S* at {y, ht{y)), so we get 



X 



-Vyht{Zt),l) 



^ and — = -Vyht{zt). 



X 



\X\ \{-Vyht{Zt),l)\ 

Making change of variables ^ = {\y, Xhtiy)) and using = A, we get 



h(t,x) 



A 



fH^'-y+-nh,iy)+t)^^^^Xy,\htiy))^\ - Uo(t,x,y) 



JU 

oo 



2^ 



JX{-x„Vyht{zt)-y+Xnht{y)+t) 



(4.17) 



JU 



X^a{t,Xy,\htiy)) 



d{\y) 



dydX 



A 

2^ 



' '^o{t,x,y)x{t,y) dydX 



JU 



x^ 



-n+l+l-l 



i^o{t,x,y)x{t,y) dydX, 



where Ko{t, x,y) = k [t ^x + V^'d{t] y, ht(y))j , and 

a(t, Xn, Xy, Xht{y)) = (a;;;^A)'a (t, x'^Ay, x~^Xht{y)) , 
and where we made a change A = x~^X in the last equality. Here also we used 



^^^^^ = A" ^x(^?l/)) where xi^iU) ^-^id all of its derivatives with respect to y are 
30unded because of (14.151) . 
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If we choose r in the definition of the cut-off function k sufficiently small, then on 
its support we have |x| ~ ~ t, and we can estimate 



(4.18) 



|Ji(t,x)| < 



J{\Zt 



X 



dX 



with 

J{X,zt) 

We will show that 
(4.19) 



Ji2'X,Zt)^\ 

,xi-Vyk.i..yy+h,iy)+t.-\^^^ x„, Xy, Xht{y))^o{t, x, y)x{t, y) dy. 



u 



|J(A,Zi)| <C(1 + A)- 



A > 0. 



Then, if we take I = n- —, and use fjiTTSl) and fjiTTOl) . we get 



\Ii{t,x)\ < Ct- — 2^ 



n — 1 n — 1 

7 



(4.20) 



(2^A)- 



A 



^ ^'l ^ 1 A> 



A-^^f - I dX = Ct 



n-l 

7 



rfA 

A-^^(A) rfA 



which is the desired estimate for Ii{t,x). 

Let us now prove fl4.19p . It will, in turn, follow from Theorem 14. 51 below. First of all 
we note that since we assumed that \d^a(t,C,)\ < Ca{C,)~'^^ for all |a| < [(n — l)/7] + l 



and A^i = n - ^ 
(4.21) 



n — 1 



+ 1, we get that 



\dya\ < C for all |a| < [{n - 1)/^] + 1. 



To write J(A, Zt) in a suitable form, we change to polar coordinates (p, u) with 
y = puj + Zt, so that 



J{X,zt) 



(4.22) 
with 

(4.23) F(p, Zt, uj) = ht{puj + Zt) - ht{zt) - pVyht{zt) ■ u, 

(4.24) f3{p, Zt, io) =a {t, Xn, X{puj + zt), Xht{puj + zt)) K^it, x, pu + Zt)x{t, puJ + zt), 

where we can assume in addition that x = unless puj + zt G U , so both p and ui 
vary over bounded sets. 

Now we can apply the following result, which has appeared in pOl for more general 
complex valued phases thus including the real-valued case of the phase function 
F in fl4.23p . The estimate fl4.19p follows from the following theorem with N = n — 1. 

Theorem 4.5 ([20]). Consider the oscillatory integral 

J(A,z/)=/ e^^*(^'^)a(x, z/)x(x) rfx , 
where N > 1, and u is a parameter. Let •j > 2 be an integer. Assume that 
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(Al) there exists a sufficiently small 5 > such that x ^ C'o°(-^<5/2(0))' where 

-85/2(0) is the hall with radius 5/2 around 0; 
(A2) $(a;, z/) is a complex valued function such that Im$(x, z/) > for all x G 

supp X O'^d all parameters u; 
(A3) for some fixed z G supp ^he function 

F{p, u, u) := $(2; + pu, u), \uj\ = 1, 

satisfies the following conditions. Assume that for each p = {00, u), func- 
tion F{-,p,) is of class C'^'^^ on supp and let us write its •y^^ order Taylor 
expansion in p at as 

7 

j=0 

where R^+i is the remainder term. Assume that we have 
(Fl) ao(/i) = ai(/i) = for all /i; 

(F2) there exists a constant C > such that Yli]=2 — ^ /^'^ /^i 

(F3) for each /i, \dpF{p,p,)\ is increasing in p for < p < 6; 
(F4) for each k < 'y + 1, dpF{p, p) is hounded uniformly in < p < 6 and p; 
(A4) for each multi-index a of length \a\ < + there exists a constant > 
such that \d^a{x, z/)| < Cq, for all x G suppx (^nd all parameters v. 
Then there exists a constant C = C^^^ > such that 

(4.25) ^)\ ^ C*(l + A) ^ for all A G [0, 00) and all parameters v. 

We refer to [20j and [22J for details. Now, the function F in fl4.23p satisfies condition 
(A3) of Theorem 14.51 because of the definition of the convex Sugimoto index 7 and 
because ht is concave. Since dyht, \a\ < 7, can be expressed via d'^^p, \a\ < 7, and 
since we have (14.161) . it also follows from (12.91) and (14. 5 p that function F satisfies 
property (F2) of Theorem 14. 5[ The proof of Proposition 14.21 is now complete. □ 

Let us now show that we can actually also insert the cut-off 1 — + \t\)C) in 
Proposition 14.21 which is necessary for the analysis of the representation (14. ip . Here 
tp G Co°°(M") such that = 1 for |^| < i, and for |^| > 1. The case of high 
frequencies |^| > 1 (for solutions) is covered by Proposition 14.21 and the proof about 
the insertion of 1 — ipii^ + 1^1)0 is similar to the proof of the following Proposition 
14. 6[ So we now restrict to t~^ < |^| < 1, since the case |^| < t^^ was covered in 
Proposition 14. 1[ 

Proposition 4.6. Let Tt, t ^ 0, he an operator defined by 

(4.26) TJ{x)= [ e^(-«+''(*^«»[l-^((l + |t|)0]a(t,0/(Orfe, 

where ^(t;^), and 7 are as in Proposition Assume that the amplitude a{t,C,) 
satisfies a{t,^) = for all \^\ > 1 and that 

|9|^a(t,0l<C„|e|-l"l for all \a\ < [{n - l)/y] + 1. 
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Let l<p<2<g<oo he such that ^ + ^ = 1. Then for t ^ we have the estimate 



(4.27) 



where N„ 



\Ttf\\Li{R'^) < C\t\ 1' 



n 



n-l 



+ 



n-l 



+ 1 



Proof. The proof of this proposition is almost the same as the proof of Proposition 
14.21 with several differences that we will point out here. Again, by interpolation, it is 
sufficient to prove estimate 



\Ttfh 



< C\t\ 



n-l 
7 



with amplitude a(t,C,) satisfying 
(4.28) |9?a(t,0l<C,|er^^-'"' 



for all |a| < [{n - + 1, 



with Ni = n-^ 



n-l 



+ 1. 



Further differences concern estimates for Ii{t,x) and l2{t,x). In general, since we 
work with low frequencies |^| < 1 only, no Besov space decomposition is necessary, 
so we do not need to introduce function \l/ and $j, so we can take \1/ = 1. 

Some additional complications are related to the fact that in principle derivatives 
of the amplitudes of operators Tt from (14.261) may introduce an additional growth 
with respect to t. In the estimate for l2{t,x) we performed integration by parts with 
operator P. Now after integration by parts the amplitude of this integral in (14.141) is 



[1 - ^((1 + mMt, 0(1 - «:) {t-'x + r Vg^(t; 0) * 



2^ 



Now, if any of the ^-derivatives falls on [1 — ^"((1 + 1^1)0]' S^t an extra factor t 
which is cancelled with t~^ in the definition of P. However, in this case we can then 
restrict to the support of Vip which is contained in the ball with radius (1 + |t|)~^, so 
we are in the situation of low frequencies |0 < t~^ again. Consequently, we can apply 
Proposition 14 . 1 1 to this integral to actually get a better decay rate of Proposition 14. 1[ 
If none of the derivatives in {P*y fall on [1 — ^/'((l + 1^1)0]' ^^e argument is the same 
as in the proof of the estimate for l2{t,x) in Proposition 14.21 

The other main difference with the proof of Proposition 14.21 is in the estimate for 
Ii{t,x). Recall now that in formula (14.171) we made a change of variables A = x~^X. 
As it was then pointed out, if r in the definition of the cut-off function k is chosen 



sufficiently small, on its support we have \xr. 



\t\. On the other hand, we have 



10 ~ A by the definition of A, since we assume that the limiting phase function cp 
is strictly positive. It then follows that (1 + |t|)^ ~ A|a;„| ~ A, and so the change 

of variables A = x~^X changes [1 — ^'((1 + 1^1)0] i^^^o " i'W] ^^e amplitude of 
/i(t, x). Justifying this argument, we can then continue as in the proof of Proposition 
14. 2[ The crucial condition for the use of Theorem l4.5l is the boundedness of derivatives 
of a in (14.211) . Here, every differentiation of a with respect to y introduces a factor 
x~^X which is then cancelled in view of assumption (I4.28p . It follows that [(n— l)/7]-|-l 
^/-derivatives of a are bounded, implying the conclusion of Theorem 14.51 This yields 
estimate (I4.27p in the way that is similar to the proof of Proposition 14.21 □ 
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Let us now turn to prove Proposition I4.4[ 

Proof of Proposition 14.41 Let us show how the proof of Proposition 14.41 differs from 

i_ 

the proof of Proposition 14. 2[ We need to prove that \I{t,x)\ < Ct , t > 0, for 
I(t,x) as in (14.131) . We note that 70 > 1, so to prove the estimate for l2(t,x) we can 
show that 1/2(^5 x)\ < Ct~^ . This can be done by integrating by parts with the same 
operator P and using (14. 101) instead of (14. 6p . As for the proof of the estimate for 
x), we can reason in the same way as in Proposition 14.21 to arrive at the estimate 
dHHD, i.e., 

with the same operator 
Ju 

Now, instead of (14.191) we will show that 

(4.29) \J{X,Zt)\ <CX'^, A>0. 

i_ 

Then, taking / = n — we get the estimate a;)| < Ct 10 in the same way as in 
estimate (I4.20p . Now, estimate (I4.29P follows from Theorem 14. 51 with N = 1. Indeed, 
let us write J(A, Zt) in the form ( ]4.22p -( l424|) with phase 

F{p, Zt, u) = htipuj + Zt) - ht{zt) - pV^ht{zt) ■ u. 

Now, by rotation we may assume that in some direction, say ei = (1, 0, . . . , 0), we 
have by definition of the index 70 that 

70 = min {fc G N : dtF{p, Zt,uj)\^,=^ ^ O}. 



Then by taking = 1 and y = cui in Theorem 14. 5[ we get the required estimate 
g29D. □ 

Now we will state the corollary of the proof of Proposition 14.41 which is similar to 
Proposition 14.61 to ensure its application to our Cauchy problem. The proof is similar 
to the proof of Proposition 14.61 

Proposition 4.7. Let Tt, t 0, be an operator defined by 

Ttf{x)= I e^(-«+''(*^«»[l-^((l + |t|)0]a(t,07(0rfe, 

where ^{t;C,) and 70 are as in Proposition 4-4' Assume that the amplitude a{t,^) 
satisfies a{t,^) = for all \C, \ > 1 and that 

\d^a{t,0\ < Cal^r'"' for all \a\ < 1. 
Let l<p<2<q<oo be such that ^ + ^ = 1. Then for t ^ we have the estimate 



(4.30) \\Ttf\\LHRn^<C\t\ -0 



where Np = (n- — + 1) (--- 



'-(---) 
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Let US finally estimate our Fourier multiplier for small t. For the analysis of very 
small t we will use the following Littlewood-Paley type theorem. 

Lemma 4.8 ([9](Theorem 1.11)). Let h = /i(^) be a tempered distribution on M" 
{n > 1) such that 

sup /''meas{^ ■ \h{^)\ > 1} < +oo, 

0</<+oo 

for some 1 < b < +oo. Then the convolution operator with J-'^^[h] is L^-L"^ bounded 
provided that l<p<2<q< +oo and ^ — = i, ^-c we have the estimate 

\\j^~~^[h]*u\\^^ <C\\u\\lp for all u G L*'(M"). 
Using this fact, we obtain: 
Proposition 4.9. Let Tt be an operator defined by 



Ttfix) 



where '&{t;C,) and a{t,C,) are as in Propositions \4-S\ or 4-4' Let n > 1, I < p < 2 < 
q < +00 and - + - = 1. Then for small t we have the estimate 



(4.31) 



where N„ = n 



\Ttf\\Li(R") < C 



In fact, the proof will yield the Besov norm B^'^ on the right hand side of the 
estimate, which is a known improvement for this type of estimates. 

Proof. In the following argument we need not use the Van der Corpt lemma, and the 
proof relies only on the Littlewood-Paley type theorem. We put K(t) = (1 + |t|)~^. 

It suffices to prove (14.311) for p ^ q, since the case p = q = 2 follows from the 
Plancherel theorem. Noting that ■(9(t; ^) and is homogeneous of order one, and making 
change of variable rj = and y = K{t)x, we get 



(4.32) ||T,/|U,(M.) = ir(t)"-? J'-' K(r^)] * J^-' [\vf^f{K{t)v) 
where we set 

mt{r]) = e'''^'^''^''''^a{t,K{t)r]){l-^{r]))\r]\'^^, 

Since a{t,^) is bounded, we have 

meas{r/ : \mt{r])\ > 1} < meas |r/ : |r/| < C^/^H-^/^A = C'^'^H-'"'^^ 



1 1 
n [ 



for each / > 0. Hence it follows from Lemma [4.81 that the convolution operator with 
rut is L^-L^ bounded, which implies that 



L'J(R") 



< c 



r 



-1 



CKity 



Ir^f^fiKm 



LP(IR") 
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where we performed the transformations K{t)r] = ^ and 



Kit) 



z in the last step. 



Thus, combining this estimate with (I4.32p . we obtain the desired estimate (I4.3ip . 
The proof of Proposition 14.91 is complete. □ 

Proof of Theorem W.W The proof of Theorem 11.11 now follows from Proposition 14 . 1 1 for 
low frequencies |^| < from Propositions 14.21 and 14.41 for large frequencies |^| > 1, 
and from Propositions 14.61 and 14.71 for intermediate frequencies < |^| < 1. We 
also use Proposition 14.91 for small times. We can note that all these propositions give 
different Sobolev orders on the regularity of the Cauchy data. 

Indeed, using representation formula for the solution established in Theorem 13.11 
we can write the solution as 

m—l 

u{t,x) = u^jt, x), 



k=0 



with 



Now, we decompose 

X) = u\{t, x) + M2(t, X) + ul{t, x) 



m 



\x 



+ E [K±(^) + 4,±(^; 0) (1 - ^((1 + mm - xiim'^^'^^^io 

with V', X e Co"^(M") such that ^(0 = x(0 = 1 for 1^1 < i and for |^| > 1. Assume 
conditions of part (i) of Theorem 11.11 Then we have estimates 



(4.33) 



\u\{tr)\\L.i^^)<C{l + \t\)-<^-^)\\f,\ 



by Proposition 14.11 
(4.34) 



with Mp = {n — {j, — ^ by Proposition 14.61 and 



(4.35) 



n-l \ 1 ' 



with = n - — + 



n-l 



the estimate 
(4.36) 



+ 1^ (p ~ g) ' Proposition 14. 2[ For small t we have 



Kit. OIU^CM") + \\ul{t, OIIl.cm") < 
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with Np = n — , hj Proposition 14.91 Putting all these estimates fl4.33p - fl4.36p 
together with similar estimates for derivatives, implies the statement of Theorem 

o □ 



As a corollary of this proof and all the propositions above, we have the following 
refinement of Theorem II. providing quantitatively different estimates for different 
frequency regions. 

Theorem 4.10. Assume fOD-f fT3|) . Let x e C'^(K) be such that xip) = 1 for 
\p\ < |j O'lT'd for IpI > 1. Let us denote 

ui = + \t\)\D\)u, U2 = {1- + \t\)\D\)xm)u, and 



u, = {l-x{{l + \t\)\D\){l-xm))u. 

Then the solution u{t,x) of fll.ip satisfies the following estimates: 
(i) Suppose that the sets 

= G R" : = 1} 

are convex for all k = 1, . . . ,m, and set 7 = max 7(S ±). In addition, suppose that 

k=l,...,m 

(1 + |t|)''<j G Li(R) forl<r< [{n - l)/7] + 1, and for all uj with \v\ + j = m. 
Let l<p<2<q< +00 and - + - = 1. Then we have the estimates 

y — — ^ p ' q 

m—l 

'1 1 ' 



^ ^ l-\-\<y.\—k^ ' 

k=0 

m—l 

\D\Dy2{t, OIU.CM") < c|tr"^(^^) IIMI^P > 1), 

' ^ i + 1 Q I + A-ip — fc ^ ^ 

k=0 
m—l 



k=0 

m—l 

|Z}^D>2(t,-)||L.(M") + P^/^>3(t,-)IU.(M") < IIMIi. _ {\t\ < 1] 

^ ^ l+la\+Np-k^ ' 

k=0 



withMp=(n-^) (i-i),iV,= (n-H-i + 

/ = 0, . . . , m — 1, anc? a an?/ multi-index. 

(ii) Suppose thafE ± is non-convex for some k = 1, . . . ,m, and set 'Jq = max 7o(S ±). 

k=l,...,m 

In addition, suppose that (1 + G L"^(M) for all iy,j with |z/| + j = m. Let 
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l<p<2<q< +00 and ^ + ^ = 1 . Then we have the estimates 

m—l 

mD^u.it, .)|U.(Rn) < C(l + J2 ll/felkf,, , jMn), (t e R), 

k=0 

m—l 

WDlDy^it, OIUnM") < C\t\Mh-.) ^ ||//clkf^^^,^^^_^(«.), > 1), 

fc=0 ^ 
fe=0 

m—l 

||L'^L>>2(i, OIIl^cm") + ||i^ti^>3(i, OIU^M") < G E ll/'^lli^ ~ (M-)' (1*1 < 1)' 

' i+|c«|+iVi,-fc^ ^ 

fc=0 

0, . . . , m — 1; anc? a any multi-index. 
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